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Eye-specific thalamic inputs converge in the primary visual cortex (V1) and form the 
basis of binocular vision. For normal binocular perceptions, such as depth and stereopsis, 
binocularly matched orientation preference between the two eyes is required. A critical 
period of binocular matching of orientation preference in mice during normal development 
is reported in literature. Using a reaction diffusion model we present the development 
of RF and orientation selectivity in mouse V1 and investigate the binocular orientation 
preference matching during the critical period. At the onset of the critical period the 
preferred orientations of the modeled cells are mostly mismatched in the two eyes and 
the mismatch decreases and reaches levels reported in juvenile mouse by the end of 
the critical period. At the end of critical period 39% of cells in binocular zone in our 
model cortex is orientation selective. In literature around 40% cortical cells are reported 
as orientation selective in mouse V1. The starting and the closing time for critical period 
determine the orientation preference alignment between the two eyes and orientation 
tuning in cortical cells. The absence of near neighbor interaction among cortical cells 
during the development of thalamo-cortical wiring causes a salt and pepper organization 
in the orientation preference map in mice. It also results in much lower % of orientation 
selective cells in mice as compared to ferrets and cats having organized orientation maps 
with pinwheels. 

Keywords: mouse VI, critical period for orientation matching, receptive field alignment, orientation selectivity, 
orientation map 



INTRODUCTION 

In the visual system, signals from the left and the right eyes first 
converge in the primary visual cortex, VI. Precise and selective 
connections between neurons are established during develop- 
mental processes for proper functioning of the nervous system. 
In simple cells in VI, the spatial layout of ON and OFF subre- 
gions determines the orientation selectivity of simple cells (Hubel 
and Wiesel, 1962). For normal binocular perceptions, such as 
depth and stereopsis, binocularly matched orientation prefer- 
ence between the two eyes is required. Little is known about 
the developmental process for binocularly matched orientation 
preference between the two eyes. Though cortical cells mostly 
have binocularly matched orientation (OR) preference between 
the two eyes in cats and monkeys, cells with interocular differ- 
ence in preferred OR (IDPO) are also reported. Blakemore et al. 
(1972) have reported a range a of ±15° (S = 6-9°) IDPOs in cat. 
Bridge et al. (2001) have reported a range of ±20° (S = 9.22°) 
IDPOs in macaque. The binocularly matched orientation pref- 
erence is established in mice after eye opening. Thus, study on 
mice provides a unique opportunity to investigate development 
of binocularly matched orientation preference between the two 
eyes (Wang et al., 2010; Sarnaik et al., 2014) during normal 
development. 

Neurons in mouse VI are selective for stimulus orientation 
(Drager, 1975; Metin et al., 1988; Van Hooser, 2007; Niell and 



Stryker, 2008). The development of orientation selectivity in 
VI has been extensively studied in cats, ferrets, and monkeys 
(Hubel and Wiesel, 1963, 1968; Chapman and Stryker, 1993; 
Sato et al, 1996; Ringach et al, 1997; White et al., 2001) in 
contrast to mice. Unlike mammals with frontally placed eyes, 
mice have more laterally placed eyes. Also in the binocular 
visual cortical area the retinal inputs from contralateral and ipsi- 
lateral eyes are not roughly equal as in cats but the ratio of 
contralateral-to-ipsilateral responses in binocular VI is approx- 
imately 2:1 (Coleman et al, 2009). Experimental studies (Drager, 
1975; Mangini and Pearlman, 1980; Metin et al, 1988; Frenkel 
et al., 2006; Niell and Stryker, 2008) suggest that in mice the visual 
cortical neurons respond to orientation stimulus similarly to that 
in cat (Tan et al, 201 1) and monkey (Van den Bergh et al, 2010). 

There is a scarcity of models on development of orientation 
selectivity in mice. Recently two computational models (Hansel 
and van Vreeswijk, 2012; Roy et al, 2013) were proposed on 
orientation selective response in layer 2/3 neurons. In both the 
models orientation selectivity was studied in layer 2/3 with feed- 
forward input from layer IV neurons and lateral connections 
within layer 2/3. Layer IV neurons are modeled with Gabor filters 
and have salt-and-pepper organization of orientation selectivity. 
At present no model exists for receptive field development in 
mouse VI, nor does a model exist that captures a salt-and-pepper 
organization of orientation selectivity. 
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Since the seminal work of Turing (1952) reaction diffusion 
equations have been applied to biological pattern formations 
extensively (Murray, 1991). Reaction-diffusion equation com- 
prises a reaction term and a diffusion term. In this paper we 
extend and apply our earlier reaction-diffusion based model 
(Bhaumik and Mathur, 2003; Siddiqui and Bhaumik, 2011) 
to model responses in mouse VI. The reaction term in our 
model captures competition for resources available at (i) the 
pre-synaptic cell, and (ii) the post-synaptic cell. The diffusion 
cooperation between neighboring cells are modeled through the 
diffusion term. We report the following. 

(i) RF development of cortical cells in monocular and binocular 
region in mouse VI: In our model cortex 39% cells are orien- 
tation selective in binocular region and 38% in monocular 
region. Reported median orientation selectivity index (OSI) 
value in mouse is 0.31 (Tan et al., 2011). For our modeled 
cells with OSI > 0.3, the mean OSIs in the two regions 
are 0.46 and 0.41, respectively. The corresponding orienta- 
tion tuning mean half widths at half height (HWHH) in 
binocular and monocular regions are respectively, 36° and 
39.4°. The mean spatial frequency in binocular zone is 0.038 
cycles/ 0 for cells with OSI > 0.3. 

(ii) Orientation (OR) preference and ocular dominance (OD) 
maps: In both monocular and binocular region we find lack 
of clustering of similar orientation preferences i.e., salt-and- 
pepper organization (Drager, 1975; Mangini and Pearlman, 
1980; Metin et al, 1988; Bonin et al, 2011) in OR preference 
map. In the binocular VI, neurons are usually contralateral 
eye dominated in mice (Drager, 1975, 1978). Wang et al. 
(2010) has reported contralateral bias in OD with a mean 
OD of 0.19 ± 0.03. The average OD in our model cortex 
is —0.1031. The negative average OD indicates that cells in 
the model binocular VI have contralateral eye bias. The OD 
map is unstructured. 



(iii) The critical period: We show how the starting and the closing 
time of the critical period affects orientation tuning, subfield 
matching in RFs, and alignment of orientation preference in 
cortical cells between the two eyes in the binocular region in 
mouse VI. 

METHODS 

THREE LAYER VISUAL PATHWAY MODEL 

Mouse VI has two zones: (i) the monocular zone, where neu- 
rons receive inputs only from the contralateral eye, and (ii) the 
binocular zone, where neurons receive inputs from both ipsi- and 
contralateral eyes as shown in Figure 1A. To obtain responses of 
cortical cells in our model mouse VI we have used a three-layer 
visual pathway model as depicted in Figure IB. In mouse the den- 
drites of retinal ganglion cells (RGCs) in the inner plexiform layer 
(IPL) of retina are separated into ON or OFF sublamina (Tian, 
2004). We have modeled ON of OFF RGC as two separate lay- 
ers. Retina for left (contralateral) eye is modeled with two 2D 
80 x 80 sheet of ON and OFF center ganglion cells lying one over 
the other. The right (ipsilateral) retina is modeled with two 2D 
40 x 80 layer of ON and OFF center ganglion cells. Mouse RGCs 
have center-surround receptive field structure with center fields 
having a radius of 5.5° and the surround field radius is 16.98° as 
reported in Grubb and Thompson (2003). Center-to-center spac- 
ing between the cells is 52' of the visual angle. The ganglion cell 
model used earlier (Wehmeier et al., 1989; Worgotter and Koch, 
1991; Somers et al., 1995; Bhaumik and Mathur, 2003) for cats is 
modified to produce realistic temporal response to visual stimuli 
in mice. The details are given in the Supplementary Material. 

The second layer models LGN. The left and the right eye 
specific LGN layers are made up of 2D 80 x 80 size and 40 x 
80 size sheets of LGN cells, respectively. In mouse dLGN ON- 
and OFF-center cells appear intermingled (Grubb et al., 2003) 
and unlike in retina do not form sublamina. It is reasonable to 
assume that initially a cortical simple cell receives inputs from 
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FIGURE 1 | The three layered visual pathway model. We have modeled 
40° of binocular and left monocular visual fields each as shown in (A). 
The left and the right eye retinae are 80 x 80 each. The left eye specific 
LGN layer is 80 x 80 with 40 x 80 region getting input from left 
monocular field and the rest 40 x 80 region from binocular field of vision. 
The right eye specific LGN layer is 40 x 80 and receives input from 
binocular field of vision. The cortical layer in the right hemisphere has 



left monocular and binocular regions of size 67 x 134 each. (B) The LGN 
layer for each eye has two sheets of cells each with center-surround 
structure — one for ON center and another for OFF center type cells. 
Each cortical cell in the left monocular region gets input from a 13x13 
section of cells from the left eye specific LGN layer, whereas each cell 
in binocular region gets input from a 13 x 13 section of cells from both 
the left eye and the right eye specific LGN layers. 
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equal number of ON- and OFF-type LGN cells. We have shown 
ON- and OFF- as two layers for ease of representation. It is 
reported that in cats each LGN cell receives strong input from 1 
to 3 retinal cells (Chen and Regehr, 2000; Jaubert-Miazza et al., 
2005). In mice about three contralateral RGCs connect to one 
LGN cell whereas one ipsilateral RGC connects to one LGN cell 
(Coleman et al., 2009). In our model we have assumed that each 
LGN cell receive input from one retinal cell. The normalization 
constant in Worgotter and Koch's (1991) model was chosen such 
that LGN cells firing rate matched experimental values (Grubb 
and Thompson, 2004) for a 90% contrast sinusoidal grating input 
to retina. The firing rate for LGN cell at the 90% contrast is 
36 spikes/s. 

The third layer models cortical layer IV of mouse VI. Third 
layer is comprised of a 2D 67 x 134 size sheet of left monocular 
cortical cells and a 2D 67 x 134 size sheet of binocular cortical 
cells. In the monocular zone each cortical cell receives synaptic 
connections from 13 x 13 left eye specific ON/OFF LGN regions 
centered at its retinotopic position. In the binocular zone each 
cortical cell receives synaptic connections from 13 x 13 left and 
right eye specific ON/OFF LGN regions. Thalamic projection 
from 13 x 13 LGN cells defines the RF of a cortical cell and cor- 
responds to inputs from approximately 21° in visual space. The 
median RF size is 21.3° in the visual space for simple cells in mice 
(Van den Bergh et al., 2010). In the model cortex, binocular zone 
covers 40° binocular visual space and monocular zone covers 40° 
left monocular visual space as shown in Figure 1A. 

CORTICAL CELL MODEL 

Membrane potential of the cortical cell is evaluated using the 
equation (Gerstner, 1998; Bhaumik and Mathur, 2003). 

ui (t) = ti (t - t fc ) +Pi J2 W IJ J2 6(t " tfL) + Rp 

J,k fL 

Where, ui (t) is the membrane potential of the cortical cell at loca- 
tion I at time t. tf c is the time when the cell produced the most 
recent spike. W^ is the strength of synapse between LGN cell at 
location J from cortical cell at location I. k is either left (1) or right 
(r). If the weight is positive, input spike train for that location is 
read from ON center LGN, otherwise from OFF center LGN. e(t - 
tfiJ is the EPSP (excitatory post-synaptic potential) generated in 
the cortical cell when there is a spike in a connected LGN cell at 
time tfL. The post-synaptic potential is scaled by the weight of the 
connection between LGN and cortical cell. The spikes occurring 
in past 90 ms are only summed up. Similarly, r|(t - tf c ) function 
is the refractory voltage of the cortical cell after it fires a spike at 
time tf c . R p is the resting potential of the cortical cell (—70 mV). 
Whenever, ui{t) exceeds threshold, a spike is generated in the 
cortical cell and the output is 1, otherwise it is 0. 

-% 

The time constants used for e 



: (t - tfl) = Fe 



g T decay 



are 5.65 ms for rise time, 2.96 ms for decay time and F = 1.35 in 
mouse compared to 12 ms rise time, 5.95 ms decay time and F = 
1.25 in cat. Also, the time constant associated with r| is 9.69 ms in 
mouse compared to 30 ms in cat. The time constants for mouse 



are smaller as compared to cat as the optimal temporal frequency 
for mouse is higher (4 Hz) compared to cat (0.4 Hz). 

We have used a modified (Bhaumik and Mathur, 2003) Spike 
Response Model (SRM) (Gerstner, 1998) to incorporate synaptic 
scaling factor /j. fi is calculated as (Bhaumik and Mathur, 2003) 

P = {H-L)f(A v )+L 

Where, / (A v ) = 1 - tanh(G (A v - 6) ) and A v = 
j ^jWjjXj (t) dt. H and L are constants and determine 
the range of the scaling factor. G and 6 are constants. A v is the 
average input the cortical cell I receives over T. Xj(t) belongs to 
{0, 1} where, Xj (t) represents the spike response of LGN cell 
at location / at time t. So, A v is related to average input spike 
rate and weights. For cells in monocular zone: H = 8.8, L = 3, 
G = 0.0061, and 9 = 180 and Binocular cells: H = 10, L = 4, 
G = 0.0061, and 6 = 100.7. Different values are used for cells 
in monocular and binocular zone to take care of difference in 
membrane potential distribution in cells in the two zones. 

For cortical cells in binocular region, the synaptic scaling 
function, fi, is calculated by taking average binocular input — 
the normal conditions of activation for such cells. f5 accounts 
for synaptic scaling mechanism, which acts over a relatively long 
time scale (over a period of several hours to days) (Turrigiano, 
2008). Therefore, for response to monocular activation of the 
mouse retina we use f3 value determined from binocular response 
of the cell. But when f$ determined from the binocular activa- 
tion was used for monocular activation, it resulted in a low value 
of the membrane potential, and consequently, very low spiking. 
But the experimental evidence suggests that monocular activation 
leads to spiking which is relatively lower compared to binocular 
one, but not too low. Recently, it is reported that inputs from 
the left and the right eye add sub-linearly in visual cortical cells 
(Longordo et al., 2013). Instead of a fixed threshold voltage, Th, 
for spiking in the cortical cell we employed a variable thresh- 
old mechanism. The variable threshold mechanism acts over the 
course of a single experiment, taking into account the varia- 
tions in input activation over successive experiments. The variable 
threshold is calculated using the following expression. 



Th = (A - B) 



De (y(A v -C)) _ e (-Y(A T -C)) 
De (Y(Av-C))_|_ e (~y(A v -C)) 



+B 



where, A = -44.5, B = 59.5, C = 100, D = 1.25; y = 0.0081, and 
the expression for A v is given earlier for f3 calculation. 

We show threshold voltage histograms in Figure 2A for cor- 
tical cells in monocular and binocular zones in the model 
cortex. The threshold histogram for binocular response peaks 
around —48 mV. This agrees with the experimentally reported 
threshold values (Tan et al, 201 1; Figure 2). The variable thresh- 
old when incorporated in the modified SRM, captures the effect 
of sub-linear addition of spike rates from monocular to binoc- 
ular experiments as reported in the literature (Longordo et al., 
2013). In Figure 2B, the binocular response spike rate is plotted 
as a function of spike rate from the sum of individual monocular 
responses. Cell response data used in Figure 2B is the response 
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FIGURE 2 | (A) Threshold value (mV) histograms for cortical cells in 
monocular (histogram on left) and binocular (histogram on right) regions 
are shown. The y-axis is the percentage of cells and x-axis is the 
threshold value in mV. (B) This is a scatter plot of observed binocular 
spike rate on y-axis and predicted binocular spike rate (as linear sum of 
individual left and right monocular spike rates) on x-axis for cortical cells 
in binocular region of cortex. The green line has a slope of 45° and the 



red line is the linear fit to data. The binocular spike rate is less than 
sum of individual monocular spike rates for higher spike rates. (C) Plots 
of p function of cells in monocular region (blue) and that in binocular 
region (red) (y-axis on left) with normalized input activation, p decreases 
as input activation increases. Plot for the threshold value (green) for all 
cortical cells with normalized input activation. Threshold value increases 
as input activation increases. 



of the cell at its preferred orientation. Data in Figure 2B is fit- 
ted with a linear function (y = 0.64 x + 7.8) and indicated with 
red line. The fitted linear line (shown in red) has a slope less 
than unity (shown as green line). At higher spike rates, the dif- 
ference between the predicted spike rates from sum of monocular 
responses and the observed binocular spike rate is more notice- 
able. We also obtain spike rates in the experimentally reported 
range of spiking under binocular and monocular activations. 
A plot of f5 and the threshold voltage as a function of pre-synaptic 
input activation is shown in Figure 2C. 

In the next subsection we modify our earlier thalamo- cortical 
synaptic weight development model (Bhaumik and Mathur, 
2003; Siddiqui and Bhaumik, 2011) for cats and obtain the 
connections between LGN and cortical cells in mouse VI. 
Biologically plausible competition and cooperation principles are 
used to model growth and decay of thalamo-cortical synaptic 
strengths. Both competition (reaction) and cooperation (dif- 
fusion) involves release of neurotrophic factors. Neurotrophins 



are activity-dependent (Bonhoeffer, 1996; Cellerino and Maffei, 
1996; Katz and Shatz, 1996; Lewin and Barde, 1996). 

THALAMO-CORTICAL SYNAPTIC WEIGHT DEVELOPMENT: MODEL 
ASSUMPTIONS 

The model is based on the following biologically plausible 
assumptions. 

1. Competition for resources available at (i) the pre-synaptic cell, 
and (ii) the postsynaptic cell: The number of synaptic connec- 
tions a LGN cell supports is constrained by its pre-synaptic 
resource. Therefore, a competition exists for a pre-synaptic 
resource where a pre-synaptic cell has a fixed amount of 
resource to distribute among its axonal branches. Similarly, a 
post-synaptic cortical cell supports a limited number of pre- 
synaptic connections depending on its post-synaptic resource. 
A competition exists between pre-synaptic LGN axons to get 
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connected to post-synaptic cortical cells. The LGN axons com- 
pete for neurotrophic factors released by the post-synaptic 
cells. Such fixed pre- and post-synaptic resources in the RGCs 
of gold fish (Hayes and Meyer, 1988) and optic tectum cells 
(Xiong et al., 1994) are reported in the literature. 
2. Diffusive cooperation between neighboring same type of ON- 
ON and OFF-OFF presynaptic LGN cells: Experimental studies 
have shown that synaptic enhancement is not restricted to 
be specific to synapses where synchronous pre- and post- 
synaptic stimulation occur. But is also accompanied by spread 
of potentiation into other inputs on the same post-synaptic 
cell (Cowan et al., 1998) within a distance of approximately 
70 (Jim (Engert and Bonhoeffer, 1997). 



THALAMO CORTICAL SYNAPTIC WEIGHT DEVELOPMENT 

In our model Wjj corresponds to the strength of the connection 
from left eye specific LGN cell at location J to the cortical cell 
at location I. It is understood in our model that positive values 
of Wjj represent connection strength from the ON center LGN 
cell, while negative values represent the strength from OFF center. 
A similar interpretation is given to Wjj which denote strengths 
corresponding to the right eye specific LGN cells. 

^ = (y[ ~ K[,) (yi - 4r) A R (I, J)a\w\j 

+ D L mfw> (i) 



The existence of both pre- and post-synaptic spread of poten- 
tiation are reported in hippocampal slice CA1 (Bonhoeffer 
et al., 1989; Bi and Poo, 2001) and hippocampal slice culture 
(Muller et al, 1995; Engert and Bonhoeffer, 1997). For thalamo- 
cortical synaptic weight development in cat in our previous work 
(Bhaumik and Mathur, 2003; Siddiqui and Bhaumik, 2011) the 
diffusive cooperation between neighboring cells models lateral 
spread of activation (Bi and Poo, 2001) both at presynaptic 
level (Bonhoeffer et al, 1989; Schuman and Madison, 1994; Tao 
et al, 2001a) and at post-synaptic level (Bradler and Barrioneuvo, 
1989; Muller et al, 1995; Engert and Bonhoeffer, 1997). In pre- 
synaptic spread of activation potentiation at synapses spreads to 
other synapses formed by the same presynaptic neuron to other 
post-synaptic neurons in the neighborhood (Bonhoeffer et al., 
1989; Schuman and Madison, 1994). The spread of potentia- 
tion on different post-synaptic cells was modeled in our earlier 
papers (Bhaumik and Mathur, 2003; Siddiqui and Bhaumik, 
2011) through the cortical diffusion term. The post-synaptic 
spread of potentiation occurs from same post-synaptic cell to dif- 
ferent presynaptic cell axons within approximately 70 u,m from 
the site of potentiation (Engert and Bonhoeffer, 1997). The 
post-synaptic lateral spread of potentiation is modeled in the 
present as well as in our earlier papers (Bhaumik and Mathur, 
2003; Siddiqui and Bhaumik, 2011) through the LGN diffusion 
term. 

Unlike in our earlier model for the synaptic weight develop- 
ment in cats (Bhaumik and Mathur, 2003; Siddiqui and Bhaumik, 
2011), we do not invoke diffusive cooperation between neigh- 
boring post-synaptic cells for synaptic weight development in 
mouse. For thalamo-cortical synaptic weight development in 
mouse we have incorporated only post-synaptic spread of poten- 
tiation. Presence of post-synaptic but no presynaptic spread is 
reported in retinotectal system in the developing Xenopus (Bi 
and Poo, 2001; Tao et al, 2001b). In retinotectal system in the 
developing Xenopus potentiation in one RGC pathway spreads 
to other retinal inputs on the same tectal neuron. But synapses 
on other adjacent tectal neurons are not affected (Bi and Poo, 
2001; Tao et al, 2001b). Such spread of post-synaptic potentia- 
tion is attributed to signaling within the post-synaptic cytoplasm 
(Bi and Poo, 2001; Tao et al., 2001b). Absence of diffusive coop- 
eration between cortical cells results in salt-and-pepper orga- 
nization in OR preference map and is discussed in the Result 
section. 



A similar equation holds for WL. The notation is explained below. 

j/j' = total presynaptic resources available at left eye specific 
LGN cell at location J. The available resource is assumed to be 
independent of J. 



+ e( - w\j) y - w 'pi» 2 J 

is the total pre-synaptic resource consumed at J. Here 0 (x) is a 
step function equal to 1 for x > 0 and 0 for x < 0. Thus, K[j gets 
contribution only from connections to cortical cells P if Wjj is 
positive (that is LGN cell is ON type) or other cortical cells P 
connected to OFF LGN cell at J if W\j is negative. The range of 
summation for P is over the two-dimensional array ofNxN 
cortical layer. 

Mechanisms behind competition between different branches 
of a single axon are reviewed in Kalil et al. (201 1). It is suggested 
that long distance signaling events mediate the transport of rate- 
limiting resources such as tubulin to the winning branch and 
thereby deprive other branches of the opportunity to grow (Butz 
et al, 2009; Van Ooyen, 2011). The factor (y\-K[j) enforces 
competition for resources among axonal branches in a left eye 
specific LGN cell. 

y 2 = total post-synaptic resources available at cortical cell at 
location I. Cortical resource available is assumed to be indepen- 
dent of I. 

K l 2I = total post-synaptic resources consumed by cortical cell 
at location I. For connection update in monocular zone in the 

cortex =(7Eq(^ q ) 2 ). 

For binocular zone, K' 2I = (JT,q (W iq ) 2 + (Wf Q ) 2 Y Here 
sum over Q is over all the M x M cells in the LGN layers. Those 
with positive values of W^q correspond to ON and those with 
negative values correspond to OFF layer. 

A post-synaptic cortical cell supports a limited number of 
pre-synaptic left and right eye specific LGN cell connections 
depending on its post-synaptic resource. The factor (y 2 — K l 2I ) 
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enforces competition among LGN cells for target space in the 
cortex. 

The arbor function, Ar(I,J) defines the region from where 
a cortical cell receives its initial unorganized thalamic afferents 
(Miller, 1994). The amount of afferents a cell receives, is deter- 
mined by the arbor window. A square window has been used for 
the results reported here. It has been shown earlier (Bhaumik and 
Mathur, 2003) that RF structure does not depend on the type of 
window used, be it trapezoidal or square. 

Aj is the activity of the left eye specific LGN cell in location "J." 
While updating a synaptic weight between a cortical cell and an 
LGN cell, the particular LGN cell has to be active, i.e., Aj = 1. If an 
LGN cell is inactive during weight update then the corresponding 
synaptic weight may decay unless helped by neighboring same- 
type cells. Before the eye opening in mouse, a dense activation 
pattern of neurons in the visual cortex is reported (Rochefort 
et al, 2009). While updating W\j we assume that cortical cell at 
location I is active. 

In Equation (1), the last term models post-synaptic lat- 
eral spread of potentiation (Bradler and Barrioneuvo, 1989; 
Muller et al., 1995; Engert and Bonhoeffer, 1997) on neighbor- 
ing LGN synapses on the same post-synaptic neuron. Engert 
and Bonhoeffer (1997) through two sets of elegant experiments, 
have shown that long-term potentiation (LTP) can spread along 
the post-synaptic dendrite from activated synapses to nearby 
synapses regardless of the history of activation of these synapses. 
The spread of LTP occurs to nearby synapses located within 
70u,m distance of activated synapse. Dijt) is the LGN diffusion 
constant and is given below. 



D L (t) 



D L , fort < t c 
0, fort>t c 



where, t c is the closing time for post-synaptic lateral spread of 
potentiation in LGN. For t > t c the synaptic weight modification 
is only governed by the reaction term. The subregion structure in 
the RF does not change — only the magnitude of weights change 
till the resources last and a steady state is reached. At 2000 epoch 
the subfields in the RF are well-formed and the RF structure does 
not depend on type of window used for Dl (t) whether it pro- 
gressively declines from Dl to zero or abruptly becomes zero at 
t = t c . 

The model parameter values in the reaction diffusion equation 
determine (i) the number of sub regions and the structure of the 
receptive filed (RF) of cortical cells and (ii) the 2D organization of 
RFs in the model cortex. The structure of RF determines the cor- 
tical cell's orientation tuning and spatial frequency selectivity. The 
2D organization of RFs determine orientation map. The effects of 
varying the model parameters, the value of Dl, the LGN resource 
y[ and the cortical resource yi are given in detail in Bhaumik and 
Mathur (2003). The number of sub fields in the RF of a corti- 
cal cell increases as Dl is reduced (see Figure 9 in Bhaumik and 
Mathur, 2003). The value of LGN resource y\ does not affect the 
structure of RF and the number of subregions (see Figure 6 in 
Bhaumik and Mathur, 2003). For low values of LGN resource y\, 
synaptic strengths are quite weak due to scarcity of resources and 
as a result the cell is not fully responsive to input stimuli. As the 



resources are increased, the synaptic weights become stronger, 
without affecting the number of sub regions and the structure of 
the RF. LGN cells compete for cortical resource y 2 . For low values 
Y2, the number of simple cells with one subregion is high. The 
synaptic strengths for cells with two or three subregions are too 
weak for < 0.5. With increase in y2>the number of cells with 
two and three subfields increases. For /2 > 1, with increase in cor- 
tical resources, the synaptic strengths increase, but the number of 
subfields remains the same (see Figure 5 in Bhaumik and Mathur, 
2003). 

We introduce a subregion correspondence factor, C ( , in 
Equation (1) for synaptic weight development for cortical cells 
in binocular zone during the critical period when subregions or 
subfields correspondence develops. During the critical period the 
evolution of connection strength Wjj is governed by the equation 



~dt~ 



(yI-K) (y2-k 1 2I )a r (i,i)a 



+ D L (t)VfW l u 



(2) 



A similar equation is used for W^. 

For updating Wy, the subfield correspondence factor C 
(Siddiqui and Bhaumik, 2011) is defined as 



-HifW^andW^ are of the same sign 
— 1 , otherwise 



(3) 



A plausible biological basis for C 1 is the N-methyl-D-asparate 
receptors (NMDARs) -mediated synaptic potentiation. A presy- 
naptic spike causes the release of glutamate, which binds to the 
NMDARs. NMDA channels allow influx of Ca ++ (Karmarkar and 
Buonomano, 2002). Larger calcium influx could produce synap- 
tic potentiation (Lisman, 1989) whereas smaller calcium influx 
could produce depression (Yuste et al, 1999). In our model, 
C 1 = +1, when at LGN location "I" synaptic connections from 
both the left and the right eyes are either ON type or OFF type. 
The active presynaptic inputs from the left and the right eye spe- 
cific LGN cells at "J," add at the post-synaptic cell and W^ grows. 
For C 1 = —1, synaptic connection from the left eye is ON (OFF) 
type but synaptic connection from the right eye is OFF (ON) 
type. Thus, both the pre-synaptic inputs are not active at the 
same time and Wjj decays unless it is supported sufficiently by 
the neighbors through the diffusive term. Had we not included C 1 
in Equation (2), irrespective of whether the synaptic connections 
from the left and the right eyes from LGN location "J" are both 
same (ON-ON or OFF-OFF) type or not, W'j grows and the left 
and the right eyes RFs of a cortical cell do not develop subregions 
or subfields correspondence i.e., the RFs in the two eyes develop 
independently. 

Synaptic plasticity in cortex is age dependent (Kirkwood et al., 
1997; de Marchena et al., 2008). Age dependent decline in visual 
cortical plasticity may be due to changes in the NMDA recep- 
tor gating properties (Carmignoto and Vicini, 1992). Further the 
mechanisms of synaptic plasticity differ across cortical layers. 



Frontiers in Systems Neuroscience 



www.frontiersin.org 



July 2014 | Volume 8 | Article 128 | 6 



Bhaumik and Shah 



Matching of binocular orientation preference 



A laminar progression of plasticity is reported in mouse visual 
cortex with critical period in layer IV lasting for a short period as 
compared to layer II/III (Jiang et al., 2007). 

The duration of NMDA-mediated excitatory post-synaptic 
currents (EPSCs) last longer in young animals as compared 
to adults. A longer lasting NMDA mediated EPSC serves to 
strengthening of synaptic connections (Carmignoto and Vicini, 
1992). It is suggested that composition of NMDARs at young 
synapses allows potentiation to occur at lower threshold (de 
Marchena et al., 2008). What factors determine the thresh- 
old for modifying synapses and how the magnitude of synap- 
tic modifications are regulated in early period of development 
in visual cortex, particularly before eye opening are mostly 
unknown. 

RESULTS 

Each cortical neuron in monocular zone receives synaptic connec- 
tions from the left eye specific LGN layers. Each cortical neuron 
in binocular zone receives synaptic connections from both the left 
and the right eye specific LGN layers. Synaptic connections are 
developed using reaction-diffusion model equations with circular 
boundary condition. The initial synaptic weights are picked from 
uniform random distribution of weights of the order of 10~ 6 . For 
the results presented in this paper the t c , the closing time for post- 
synaptic lateral spread of potentiation in LGN is 3000 epochs. 
LGN diffusion constant value Dl = 0.0055 for t < t c . The cortical 
resource, yi_ = 1.5. Contra-lateral eye makes stronger connec- 
tions as compared to ipsi-lateral eye (Coleman et al, 2009). Thus, 
in our model LGN resources of contra-lateral eye (Left eye here), 
y[ =1.2 and that of ipsi-lateral eye (Right eye here), y{ = 0-8. 

In the monocular zone, synaptic connections for cortical 
cells are developed using reaction-diffusion equation given in 
Equation (1). In absence of diffusive spread of potentiation on 
neighboring LGN synapses on the same post-synaptic neuron 
i.e., with the LGN diffusion constant Dl = 0, the synaptic weight 
development is governed by only the reaction term and developed 
RFs have random distribution of ON and OFF synaptic weights. 
No ON or OFF sub field develops in absence of LGN diffusion. 
When both reaction and diffusion term are included for devel- 
opment of synaptic weights as in Equation (1), from the initial 
random distribution of synaptic weights at f = 0, the relatively 
higher values of synaptic weights acts as nucleation centers and 
at around epoch 100, the left and the right RFs of the cortical 
cell develop small patches of ON or OFF subregions. The for- 
mations of patches occur due to cooperation among ON (OFF) 
synapses helping other neighboring ON (OFF) synapses to grow 
and push out any OFF (ON) synapses existing in a patch. The 
cooperation phenomenon is gradual and is due to diffusion in 
the LGN. At epoch 3000, RFs have well-defined segregated ON 
and OFF subregions with gradual transition from ON (OFF) sub- 
region to OFF (ON) subregions. At t = t c = 3000 epochs, the 
utilization of resources are almost saturated, and the change in 
synaptic weights are negligible. For t>3000 epochs the synaptic 
weight modification is only governed by the reaction term. The 
subregion structure in the RF does not change — only the magni- 
tude of weights change till the resources last and a steady state is 
reached. 



To achieve computational speed, we assume LGN cell and its 
neighbors are active during weight update for the results pre- 
sented in this article. In the Supplementary Material we have 
shown that the nature of RFs developed with spontaneous activity 
in LGN remains qualitatively same. 

Critical period is the period of heightened plasticity during 
which the orientation alignment between the left and the right 
eye receptive fields for cortical cells in the binocular zone takes 
place in normal development. Wang et al. (2010) reported that in 
mice, the left and the right eye receptive fields have poor match- 
ing of preferred orientation at the start of the critical period but 
they become matched by the end of it. Let t = C, ter be the epoch 
at which the critical period starts and t = t c be the epoch for the 
end of critical period for orientation plasticity. In binocular zone 
therefore, we allow the left and the right eye receptive fields to 
develop independently till 500 epochs using Equation (1). From 
500 to 3000 epochs we develop synaptic weights for binocular 
cells using Equation (2) wherein subregion correspondence fac- 
tor is included in the reaction diffusion equation. The left and the 
right eye RFs at 500 and 3000 epochs for a 5 x 5 section of cor- 
tex in the binocular zone are shown in Figure 3. According to the 
feed-forward model proposed by Hubel and Wiesel (1962), ori- 
entation selectivity arises from specific arrangement of geniculate 
inputs. Ours is a feed forward model and the preferred orienta- 
tion of modeled cortical cells is determined by the layout of the 
elongated ON and OFF subregions as shown in Figure 3B. RFs are 
locally highly diverse, with nearby neurons having largely dissim- 
ilar receptive fields. The ON and the OFF subregions are shown in 
Gray-scale with white (black) color representing strong synaptic 
connection from ON (OFF) LGN cells. The shading is propor- 
tional to the strength of the ON/OFF synaptic connections from 
the LGN cells. The choice of 500-3000 epochs as time window 
for critical period for orientation plasticity is justified later in this 
section. 

CELL RESPONSE 

To obtain cell response using SRM moving sinusoidal gratings 
at 90% contrast (Grubb and Thompson, 2003) were used. The 
velocity of gratings was fixed at 100°/s orthogonal to the orien- 
tation. This results in temporal frequency of 4Hz for sinusoidal 
grating with spatial frequency of 0.04 cycles/ 0 . Stimulus orien- 
tation was varied from 0° to 180° in steps of 18°. Peri-Stimulus 
time histogram (PSTH) was made for each response by binning 
the response in bins of 100 ms each and calculating the spike rate 
in each of these bins. Spike rates per second were computed for 
individual bins and the response was then averaged over the 30- 
recorded Peri-Stimulus time histograms. The cell spike response 
for any given orientation of input stimulus is the maximum 
response obtained in the averaged histogram. The mean spike rate 
for binocular cells is 22 spikes/s and for left monocular cells it is 
19 spikes/s. This agrees with data in literature (Metin et al, 1988). 

Approximately 60% non-orientation selective cells are 
reported in mouse VI (Mangini and Pearlman, 1980; Metin 
et al, 1988; Hiibener, 2003) i.e., around 40% cells are orientation 
selective. In our study, in the binocular zone we observed 
39% orientation selective cells when stimulated binocularly. In 
monocular zone 38% cells are orientation selective. 



Frontiers in Systems Neuroscience 



www.frontiersin.org 



July 2014 | Volume 8 | Article 128 | 7 



Bhaumik and Shah 



Matching of binocular orientation preference 





Left and Right RFs at 500 









H 










tm 








li 








□ 


















Si 




Left and Right RFs at 3000 



FIGURE 3 | (A,B) The left and the right RFs for a 5x5 section of cortex 
in the binocular zone are shown at 500 epochs, and 3000 epoch. RFs 
are locally diverse where near-by neurons have largely dissimilar receptive 
fields. The ON and OFF subregions are shown in Gray-scale with white 
(black) color representing connections from ON (OFF) LGN cells. The 
shading is proportional to strength of connection from LGN cells. At 500, 
epochs the sub-field structure is found only in few cells. The left and the 



right RFs are not similar as sub region correspondence factor is not 
included before 500th epoch. At 3000, the sub-field structure is visible. 
Due to sub region correspondence factor acting between 500 and 3000 
epochs, the left and the right RFs become similar. Black boxes are 
marked around RFs for cells which are oriented in left and right eyes 
individually and the orientation preference difference between the left and 
the right eyes is less than 30°. 



Some studies (Metin et al., 1988; Niell and Stryker, 2008) 
report that the tuning width of mouse cells is very good, around 
20° while others (Metin et al., 1988; Van Hooser, 2007) report that 
they could not get many highly tuned cells when inhibition was 
turned off. In our study we have feedforward excitatory connec- 
tions from LGN to cortex. We get some cells that are highly tuned, 
while there is also a considerable number that have large HWHH 
and hence, are poorly tuned. Cells with three subregions shows 
better tuning as compared to cells with two subregions, particu- 
larly if one of the two subregion is bigger in size. Also cells with 
patchy subregions i.e., not well-formed subregions are poorly ori- 
entation selective. For cortical cells in binocular region, HWHH 
histogram for right monocular, left monocular, and binocular 
responses are shown in Figure 4A. A number of cells have low 



values of HWHH (<20°) in their monocular response. However, 
only those cells that do not have large difference between the 
left and the right eye orientation preference have low values of 
HWHH in binocular response. HWHH histogram for cortical 
cells in monocular region is shown in Figure 4B. Note that the 
HWHH histogram for cells in monocular region is similar to the 
HWHH histogram for binocular response of cells in binocular 
region. This is because spike rates for monocular cells are sim- 
ilar to the spike rates for binocular response in binocular cells. 
Figure 4C depicts spike rate as function of HWHH for cells in 
binocular zone. We observe that for low spike rate, we are more 
likely to get better orientation tuning. 

We have also obtained OSI of the modeled cells. OSI is cal- 
culated by making a vector for each orientation response with 
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FIGURE 4 | (A) HWHH histogram of the left monocular, the right monocular, 
and the binocular responses for cortical cells in a 35 x 60 section inside the 
binocular region. There is a large chunk of cells which are poorly orientation 
tuned, but there is also a large number of cells with good orientation tuning. 



Tuning properties of monocular responses are better than binocular responses. 
(B) HWHH histogram for a 35 x 60 section of cells in the monocular region. (C) 
Scatter plot between HWHH and maximum spike rate (spikes/s) of binocular 
response for a 35 x 60 section of cells in the binocular region of the cortex. 



magnitude equal to the magnitude of response and angle being 
twice the stimulus orientation. The OSI is calculated using the 
relation (Tan et al., 2011) 

J( E e Resin(29)) 2 + ( £ e R e cos(29)) 2 
OSI= ± — 

where R§ is the response for stimulus orientation angle 9. OSI 
value varies between 0 and 1. Median OSI value in mouse is 0.31 
(Tan et al, 2011). In the model cortex there is a large chunk of 
cells, which are poorly orientation tuned (OSI < 0.3), but there 
is also a large number of cells that show good orientation tuning 
(OSI > 0.3). This result is nearer to that observed in Hagihara 
and Ohki (2013). 

Reported optimal spatial frequency range across the cell pop- 
ulation in mice ranges from 0.02-0.09 cycles/ 0 (Gao et al., 2010) 
to 0.003-0.1 cycles/ 0 (Van den Bergh et al, 2010; Vreysen et al, 
2012). We have studied the spatial frequency response of our 
modeled cells. For obtaining spatial frequency response a cen- 
trally located cortical section in binocular region was chosen. 
Further, cells (N = 374) that are orientation selective in their 
binocular response and have difference in the preferred orien- 
tation, AOR < 30°, between left eye and right response were 
selected. The spatial frequency responses of these cells were 
obtained by exposing the cells to drifting gratings at preferred 
orientation as obtained from their binocular response, at dif- 
ferent spatial frequencies but at same velocity (100°/s). The 
spatial frequency responses of three sample cells are shown in 



Figure 5A. The histogram for the optimal spatial frequency is 
shown in Figure 5B. For cells (N = 138) with at least mod- 
erate tuning (OSI > 0.3) the spatial frequency histogram is 
shown in Figure 5C. The mean preferred spatial frequency of 
these cells is 0.038 cycles/ 0 and close to the values reported 
(Gao et al., 2010; Van den Bergh et al., 2010; Vreysen et al., 
2012). 

OR PREFERENCE AND 0D MAP 

Two orientation preference maps for 32 x 32 section of cor- 
tex inside binocular region and monocular region are shown in 
Figures 6A,C. The white color lines depict the binocular preferred 
OR of cells. No orientation line is shown for the cells that are not 
orientation tuned. The orientation preferences of nearby cortical 
cells are uncorrelated and we observe a salt and pepper OR pref- 
erence map. In cat visual cortex, neurons are arranged according 
to their preferred orientation resulting in smooth OR preference 
map with pinwheels. In mice there is no evidence of such smooth 
OR preference map (Drager, 1975; Mangini and Pearlman, 1980; 
Metin et al., 1988; Bonin et al, 2011) although neurons are tuned 
to orientation of the visual stimulus. In a salt-and-pepper orga- 
nization of OR preference map, the receptive fields with different 
preferred orientations could sample a visual scene uniformly. In 
cats and ferrets with structured OR preference maps at certain 
visual field positions some degree of overrepresentation of certain 
orientations are likely. But in combination with eye movements 
this disadvantage may be overcome (Kaschube, 2014). 

Hubel and Wiesel (1963) had suggested that smooth OR 
preference map in cat visual cortex may lead to wiring length 
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FIGURE 5 | (A) The spatial frequency response of 3 sample cells. The 
response at sample frequencies (0.01-0.06 cycles/" with intervals of 0.01 
cycles/") and velocity of 100°/s is fitted with cubic spline to determine 
optimal spatial frequency. Optimal spatial frequency for Cell 1— Cell 3 are 
0.042 cycles/", 0.022 cycles/", and 0.048 cycles/", respectively. (B) The 



Optimal spatial frequency histogram for all cells (A/ = 374) from a patch 
in binocular region that are oriented in left and right RFs and have AOR 
< 30°. (C) Optimal spatial frequency histogram for ceils from (B) that 
also have moderate orientation tuning (OSI > 0.3). Mean in (C) is 0.038 
cycles/". 



minimization in cortex. More recently it is proposed that the dif- 
ferences in the structure of orientation maps in different species 
may be explained by differences in intra-cortical circuits and 
wiring length minimization (Koulakov and Chklovskii, 2001). 
One of us has shown (Bhaumik and Mathur, 2003; Siddiqui and 
Bhaumik, 2011) that local diffusive interaction among cortical 
cells results in formation of smooth OR preference map with 
pinwheels. Lack of diffusive cooperation among neighboring cor- 
tical cells results in salt and pepper organization of OR preference 
maps as shown in Figures 6A,C. We, therefore, suggest that a lack 
of columnar organization in OR in rodents is due to absence of 
diffusive cooperation among neighboring cortical cells in these 
animals during development. 

Formation of RFs with ON/OFF subregions in orientation 
selective visual cortical cells is determined by local diffusive 
interaction among neighboring LGN cells. Hence development 
of orientation selectivity takes place irrespective of presence or 
absence of columnar organization of orientation selectivity. The 
reaction diffusion model used by us in this paper for mouse and 
in our earlier work on cat (Bhaumik and Mathur, 2003; Siddiqui 
and Bhaumik, 2011) can be used to study the formation of RF 
of orientation selective cells in diverse animals such as mice or 
rats with no OR preference map to ferrets and cats with smooth 
OR preference map with pinwheels. Orientation preference his- 
togram for binocular cells and monocular cells are shown in 



Figures 6B,D. All orientations are almost equally present. 39% 
of cells in binocular zone in our model cortex are orientation 
selective. 

OD for cells in binocular region is a measure of dominance of 
response of one eye over the other. OD value for our modeled cells 
is calculated as: 

L 

OD = where, L = R l( e ) and R = Rr ^ 

Ri(9) is the response when stimulus is applied to only the left 
eye and R r (6) is the response to stimulus applied only to the 
right eye. 8 is stimulus orientation angle. OD histogram for cells 
in binocular zone is shown in Figure 6E. OD lies between —1 
and +1. Negative OD values correspond to left eye dominance 
and positive OD values correspond to right eye dominance. OD 
values near 0 indicate binocularity. The cell population shows 
contra-lateral dominance with mean OD of —0.1031. The OD 
map for the same 32 x 32 section of cortex considered earlier for 
Figure 6A, is shown in Figure 6F. The OD map is unstructured. 

CRITICAL PERIOD AND ORIENTATION PLASTICITY 

In mouse, the critical period starts at around P21 and lasts till P35 
(Huberman et al, 2008). Mice aged between P20-P23 and P31- 
P36 are used to study binocular orientation preference matching 
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FIGURE 6 | (A,C) Two orientation preference maps for a 32 x 32 section 
of cortex inside binocular region and monocular region, respectively. The 
lines depict orientation preference for cells. For cells that are not 
orientation tuned, no line is shown. We observe a salt and pepper 
orientation preference map. (B) Orientation preference histogram for 



binocular cells. (D) Orientation preference histogram for monocular cells. 
All orientations are almost equally present. (E) OD histogram of the cells 
in binocular region. Mean OD is —0.1031 and depicts contra-lateral 
dominance. (F) OD map for the section of the cortex shown in (A). The 
OD map is unstructured. 



(Wang et al., 2010). In our study 500 epoch and 3000 epoch 
corresponds to start and end of orientation plasticity. There is 
abundance of literature on ocular dominance plasticity (ODP). 
Mechanisms such as long-term depression (LTD), LTP and neu- 
rotropic signaling are proposed for ODP (Espinosa and Stryker, 
2012). Literature on Orientation plasticity is limited (Wang et al., 
2010; Sarnaik et al., 2014) and no mechanism has been proposed. 

For binocular RF development prior to the critical period, we 
use Equation (1) till 500 epochs. We start the iteration with very 
small (of the order of 10~ 6 ) random values of synaptic weights. 
Due to reaction term, the weights quickly increase and observ- 
able responses from cells with orientation tuning start from 500 
epochs onward. A significant number (approximately 20%) cells 
display orientation tuned monocular responses at 500 epochs. 



But as the receptive field in two eyes develop independently, only 
4% cells have both tuned left and right monocular responses 
with absolute value of |OD|< 0.8 and maximum spike rate in 
left and right monocular responses greater than 3 spikes/s (base- 
line firing rate of cortical cells). Experimentally oriented cells 
are observed at the start of critical period (Wang et al., 2010). 
In our simulation we take 500 epochs as the start of the criti- 
cal period. We develop synaptic weight using Equation (2) from 
500 to 3000 epochs. Histogram of HWHH of binocular cells is 
shown in Figure 7 A starting from 500 epochs to 3000 epochs at 
an interval of 500 epochs. With increase in epochs the number 
of binocular orientated cells increases from 4% at 500 epochs to 
39% at 3000 epochs. The increase in number of cells with bet- 
ter binocular tuning is due to the improvement of orientation 
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FIGURE 7 | (A) HWHH histogram for the cells in a 35 x 60 patch 
from the binocular region in the cortex. The histograms are shown 
from 500 to 3000 epochs at an interval of 500 epoch. (B) 
Histograms of preferred orientation difference between the two eyes 
at (i) 500, (ii) 1500, and (iii) 2500 epochs. In (iv) the histograms 



in (i)-(iii) are compared with uniform distribution (green line). 
Uniform distribution is the expected distribution of |AOR| when the 
left and right receptive fields develop independently. At 500 epochs, 
|AOR| distribution is almost uniform. There is not much improvement 
in |AOR| between 1500 and 2500 epochs. 



tuning in monocular responses as well as improved alignment of 
orientation preferences between the two eyes. 

During the critical period the alignment of the left and the 
right eye preferred orientation takes place. Figures 7Bi— iii depicts 
the histograms of preferred orientation difference between the 
two eyes at 500, 1500, and 2500 epochs. With increasing num- 
ber of epochs, the number of binocular orientated cells increases 
along with improvement in alignment between preferred orien- 
tations of the two eyes. Figure 7Biv shows a plot of cumulative 
distribution of the percentage of cells as a function of preferred 



orientation difference between the two eyes, AOR, from 500 
epochs to 2500 epochs at an interval of 1000 epochs. The green 
line in the plot indicates uniformly distributed orientation differ- 
ence. Note that at 500 epochs, the orientation difference distribu- 
tion is almost uniform. This is expected because till 500 epochs 
the RFs of the two eyes developed independently. The preferred 
orientation alignment improves significantly up to 1500 epochs; 
thereafter, it is slower. We take 3000 epochs as the closing of the 
time window for orientation plasticity in juvenile mouse. How 
the choice of starting and closing time for the time window for 
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orientation plasticity affects the preferred orientation alignment 
between the two eyes and the binocular orientation tuning is 
discussed next. 

START AND END OF CRITICAL PERIOD 

The choice for deciding on the starting time (C-iter) for the crit- 
ical period for orientation plasticity is guided by the following 
three factors. 

1. At the beginning of the critical period there must be some 
cells that are orientation tuned in both eyes so that the pre- 
ferred orientation for the two eyes can be ascertained. Also, 
the preferred orientation in most of the cortical cells should be 
mismatched between the two eyes. 

2. Similarity of orientation preference between the two eyes 
should improve during the critical period. At the end of crit- 
ical period (3000 epoch in our case) we must have similar 
amount of alignment of preferred orientation of the two eyes 
as reported in juvenile mice. 

3. At the end of plasticity time window approximately 40% cells 
should be orientation tuned as reported in literature. 

Histogram of preferred orientation between two eyes, A OR, at 
3000 epochs for C-iter = 0, 500, and 1500 epochs are shown in 
Figures 8Ai-iii. We have plotted in Figure 8Aiv the cumulative 
distribution of AOR, at 3000 epochs for three different C-iter val- 
ues. The green line depicts uniform distribution. Best orientation 
alignment for the two eyes are obtained for C-iter = 0 epoch. 
However, at 0th epoch no cells are orientation tuned, thus first 
condition for choice of C-iter is not satisfied. For C-iter = 1500 
epochs the alignment in preferred orientation is poor and the 
cumulative distribution is almost similar to random distribution. 
With C-iter = 500 epochs, at 3000 epoch we get 39% orientation 
selective cells in the binocular zone with orientation alignment 
similar to the ones reported in literature (Wang et al, 2010). The 
HWHH of cells for binocular input at maturity for three differ- 
ent C-iter values are shown in Figure 8B. The mean HWHHs are 
51.45°, 57.27°, and 58.41° at C-iter = 0, 500, and 1500 epochs. 
Note that the mean HWHH data contains all orientation selective 
cells i.e., cells with HWHH < 90°. If we consider only moderate to 
high orientation selective cells i.e., cells with OSI > 0.3 the mean 
HWHH improves. For instance for C-iter = 500, at 3000 epoch 
mean HWHH is 36° for cells with OSI > 0.3. 

The end of critical period for orientation plasticity is deter- 
mined by two factors. 

1. Orientation tuning, and 

2. Alignment of preferred orientation between the two eyes. 

In the result presented in this paper 3000 epochs is taken as the 
time for end of critical period. To study the effect of extending the 
duration of critical period on RFs, we carried on synaptic weight 
development for different values of t C) ranging from 3000 epoch 
to 70,000 epoch. Figure 9A shows left and right eye RF for some 
sample cells from the model cortex for differenr values of t c . For 
larger values of t c , the diffusion continues for longer time and 
the sub-field structure of the receptive fields is lost. Consequently 



orientation selectivity is also lost. The critical period of plasticity 
in layer IV in mouse is shorter than that in layer II/III (Jiang et al., 
2007). Our study shows that a critical period beyond 2500-3000 
epoch in layer IV causes loss in orientation selectivity. 

The alignment of preferred orientation comes from subregion 
correspondence in the RFs of the left and the right eyes. To quan- 
tify subregion correspondence we calculate a matching index for 
RFs for the left and the right eyes. We define the matching index 
for a cortical cell as the fraction of the weights in left and right 
receptive fields that have same sign (either ON-ON or OFF-OFF). 
At each epoch we obtain the matching index for each cortical cell. 
Next we determine the overall matching index. Overall matching 
index at an epoch is the matching index averaged across all corti- 
cal cells. Matching index is only a measure of subregion/sub-field 
correspondence and not of orientation alignment. In Figure 9B, 
we plot the change in overall matching index when RFs were 
developed with three different values C-iter = 500, 1 500, and 2500 
epochs. In all the three cases, the matching index is 0.5 before the 
critical period starts as the left and the right receptive fields for a 
cortical cell develop independently. After the critical period starts, 
the matching index increases, as expected. It takes more iterations 
for matching index to reach similar value in case of C-iter = 2500 
compared to C-iter = 500, but orientation tuning is lost as the 
number of epoch increases and we observe receptive fields with 
blob like single region. The rate of increase in matching index is 
faster when critical period starts earlier. 

Our study predicts that the starting and closing time for critical 
period during normal development affect the orientation prefer- 
ence alignment between the two eyes and orientation tuning in 
cortical cells. Therefore, animals with similar OR preference map 
organization but different developmental time scale are expected 
to have different distribution of AOR and percentage of OR tuned 
cells. This prediction needs experimental verification. 

DISCUSSION 

We have used a reaction diffusion model to develop RFs to study 
orientation selectivity in mouse visual cortex during the criti- 
cal period when alignment of orientation preference between the 
two eyes takes place. Development of orientation selectivity in 
both monocular and binocular regions are modeled for mice in 
this paper. For cats we had earlier modeled development of ori- 
entation selectivity in binocular region (Siddiqui and Bhaumik, 
2011, 2013). In the model mouse cortex binocular zone covers 
40° binocular visual space and monocular zone covers 40° left 
monocular visual space. To obtain cortical cell response using the 
three layer visual pathway model we have (i) implemented a RGC 
model for mouse, (ii) developed a SRM for cortical cell to cap- 
ture realistic spike rates and responses in binocular as well as in 
monocular region of the mouse cortex, (iii) presented thalamo- 
cortical synaptic weight development in mice, and (iv) studied the 
matching of binocular orientation preference during the critical 
period. 

Our model captures the diversity in RFs and orientation pref- 
erences (see Figure 3) in local cell population as reported in 
Bonin et al. (2011). The local diversity in orientation prefer- 
ence in mice is captured in the salt and pepper orientation maps 
(see Figures 6A,C). Approximately 75% of cells in adult ferret 
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FIGURE 8 | (Ai-iii) Histogram of preferred orientation difference between 
the two eyes, AOR, at maturity (3000th epochs) for C-iter = 0, 500, and 
1500 for cells from the same 35 x 60 patch of cortex used for Figure 7. 
(iv) The cumulative distribution of AOR, at 3000 epochs is shown for 



three different values of C-iter. The green line depicts uniform 
distribution. (B) The HWHH histograms for binocular oriented cells from 
the same patch of cells as in (A) for (i) C-iter = 0, (ii) C-iter = 500, and 
(iii) C-iter =1500. 



cortex (Chapman and Stryker, 1993) and about 90% of cells 
in adult cat cortex (Bishop and Henry, 1972) are orientation 
selective. In mouse VI around 40% cells are orientation selec- 
tive (Mangini and Pearlman, 1980; Metin et al, 1988; Hiibener, 
2003). We attribute the lower % of orientation selective cell in 
mice as compared to ferrets and cats to absence of diffusive coop- 
eration between neighboring cortical cells during development. 
Absence of diffusive cooperation between cortical cells also causes 
salt and pepper orientation map in mice. To get a quantitative 
understanding of how diffusive cooperation among near neigh- 
bor cortical cells during development affects the % of orientation 
selective cells and tuning properties we studied development of 



orientation selectivity in our model cat cortex (Siddiqui and 
Bhaumik, 2011) (i) in absence, and (ii) in presence of diffu- 
sive cooperation among cortical cells. The orientation tuning 
histogram for the two cases are shown in Figures 10A,B> respec- 
tively. All parameters for the developing RFs in cat for Figure 10B 
are same as given in our earlier work (Siddiqui and Bhaumik, 
2011). In absence of diffusive cooperation among cortical cells, 
57.12% cells are orientation selective as compared to 84.34% 
when diffusive cooperation is present. 

It was suggested that cortical neurons with IDPOs forms "sec- 
ond neural mechanism for depth perception" (Blakemore et al, 
1972). Bridge and Cumming (2001) and Bridge et al.'s (2001) 
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FIGURE 9 | (A) The effect of extending the duration of critical period on 
RFs. The left eye and the right eye specific RFs for some sample cells 
are shown. The RFs shown on the left of the figure are for the left eye 



and those shown on the right are for the right eye. (B) The change in 
overall matching index for three different values of C-iter — 500, 1500, and 
2500 epochs, respectively. 
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FIGURE 10 | HWHH histogram in model cat cortex, (A) in absence, and (B) in presence of diffusive cooperation among cortical cells. 



study did not support this hypothesis and they showed that VI 
cells with IDPO in monkey are not effective in horizontal slant 
detection. We have recently shown (Siddiqui and Bhaumik, 2013) 
that cortical cells in VI do not possess slant selectivity. During 
normal development binocularly matched orientation preference 
between the two eyes develops. Cells with IDPO are cells with 
imperfect binocular matching at the end of the critical period 
for the alignment of the left and the right eye preferred orien- 
tation. We have shown in Table 1 the distribution of cells with 
IDPO in the model cat visual cortex with and without diffusive 



cooperation among cortical cells. There is no significant effect of 
cortical diffusion on orientation preference matching. 

We have modeled development of orientation selectivity 
through connections from thalamus to cortical simple cells in 
layer IV in mice VI in the present paper as well as in our earlier 
work for layer IV in cat VI (Bhaumik and Mathur, 2003; Siddiqui 
and Bhaumik, 2011). Orientation selectivity emerges from spe- 
cific connections from thalamus to cortical simple cells in layer IV 
(Hubel and Wiesel, 1962; Reid and Alonso, 1995, 1996; Ferrester 
et al., 1996; Kara et al., 2002; Ohki and Reid, 2007). Orientation 
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Table 1 | Distribution of IDPO with and without cortical near neighbor 
interaction. 
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8.19 


30° < |IDPO| <40° 


17.34 


7.77 



selectivity proceeds, or emerges before the clustering of hor- 
izontal projections (Sur and Leamey, 2001). At P23 approxi- 
mately one quarter of ferret cortical neurons studied by Chapman 
and Stryker (1993) displayed orientation selectivity. Significant 
Clustering of horizontal connections were seen around P27 
(Ruthazer and Stryker, 1996). Recently it is reported that in mice 
emergence of orientation selectivity does not depend on the pre- 
cise arrangement of local horizontal connections (Ko et al., 2013). 
Neurons first acquire feature preference from feedforward inputs 
before the onset of sensory experience (Ko et al, 2013). At the 
eye opening neurons are already highly selective for visual stimu- 
lus. After eye opening, local connectivity reorganizes extensively. 
The nascent orientation map at the eye opening is likely to serve 
as substrate for the reorganization. Local horizontal connections 
in both in animals having organized orientation preference map 
as in cats and salt-and-pepper map as in mice are approximately 
isotropic (Burkhalter and Charles, 1990; Bosking et al, 1997; 
Buzas et al., 2006). Moreover, local connections in both cats 
(Martin and Schroder, 2013) and mice (Ko et al, 2011, 2013) 
more frequently link neurons with similar preferred orientations. 

The mouse has emerged as an important model system for 
studying development of neural circuitry. This paper, to the best 
of our knowledge, presents for the first time a model on the devel- 
opment of receptive field for orientation selective cortical cells in 
layer IV in mouse VI. The model captures lack of clustering of 
similar (i) orientation selectivity preferences and (ii) OD in mice. 
Our study shows how the starting and closing time for critical 
period affect the orientation preference alignment between the 
two eyes and orientation tuning in cortical cells. The reaction 
diffusion model can be used to capture the formation of RF of 
orientation selective cells and their spatial organization in diverse 
animals such as mice or rats with no organized OR preference 
map to ferrets and cats with smooth OR preference map with 
pinwheels. 
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Figure 1S | (A) Synaptic weight development according to the LGN activity 
pattern when LGN cell at location J and its neighbors are active during 
weight update. (B) Synaptic weight development when LGN layer is 
activated by spontaneous activity pattern. 
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